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It is clear that optical selection effects have distorted the "true" GRB redshift distribution to its 
presently observed biased distribution. We constrain a statistically optimal model that implies 
GRB host galaxy dust extinction could account for up to 40% of missing optical afterglows and 
redshifts in z = — 3, but the bias is negligible at very high-z. The limiting sensitivity of the 
telescopes, and the time to acquire spectroscopic/photometric redshifts, are significant sources of 
bias for the very high-z sample. We caution on constraining star formation rate and luminosity 
evolution using the GRB redshift distribution without accounting for these selection effects. 



Gamma-Ray Bursts 2012 Conference -GRB2012, 
May 07-11, 2012 
Munich, Germany 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Lieence. 



http://pos.sissa.it/ 



The mystery of the missing GRB redshifts 



David M. Coward 



1. GRB optical selection effects 

We define gamma ray burst (GRB) optical afterglow (OA) selection effects as the combination 
of sensitivity limited optical follow-up and phenomena, astrophysical and instrumental, that reduce 
the detection probability of an OA. Some of the more widely understood effects are discussed by 
Fynbo et al.[jl|] in detail. Optical biases have reduced the fraction of Swift triggered OAs, and have 
introduced a selection towards detecting the brightest OAs, hence the more nearby bursts. Addi- 
tionally, there are biases that distort the redshift distribution over certain redshift ranges (see e.g. 

IBB): 

1. Malmquist bias: This bias arises because the telescopes and instruments acquiring OA ab- 
sorption spectra (and photometry) are limited by sensitivity. In reality, the instruments ac- 
quiring redshifts are biased to sampling the bright end of the OA luminosity function. To 
account for this bias, it is necessary to have some knowledge of OA luminosity function 
(which is uncertain especially at the faint end), and an estimate of the average sensitivity 
limit of the instruments. This is the most fundamental bias that encompasses all flux limited 
detection and is the basis for modelling a selection function for OA/redshift measurement. 

2. Redshift desert: The so-called redshift desert is a region in redshift (1.4 < z < 2.5) where it 
is difficult to measure absorption and emission spectra. As redshift increases beyond z ~ 1, 
the main spectral features become harder to recognize as they enter a wavelength region 
where the sensitivity of CCDs starts to drop and sky brightness increases. Beyond z ~ 1 .4, the 
spectral features move beyond 1 /im, i.e., into the near-IR. In the case of actively starforming 
galaxies at z > 1.4, these are several narrow absorption lines over the UV continuum, most 
of which originate in the ISM of these galaxies. 

3. Different redshift measurement techniques: Historically, because of the deficiency in pre- 
Swift ground-based follow-up of GRBs, there was a strong bias for imaging the brightest 
bursts. Because the brightest bursts are predominantly nearby, a significant fraction of the 
first GRB redshifts were obtained by emission spectroscopy of the host galaxy. In the Swift 
era (from 2005 onwards), an optical afterglow (OA) is usually required to measure a redshift. 
For most high-z GRBs, this is achieved by absorption spectroscopy of the GRB afterglow. 
The host galaxies are usually too faint to make a significant contribution to the spectra. Most 
GRB spectroscopic redshifts are acquired by large aperture ground based telescopes, includ- 
ing VLT, Gemini-S-N, Keck and Lick (see [[|] for a more complete list along with specific 
spectroscopy instruments). The measurement of a GRB redshift depends strongly on the 
limiting sensitivity and spectral coverage of the spectroscopic system. This bias is expected 
to manifest at high-z, where the optically brightest OAs are near the limiting sensitivity of 
the telescope. 

4. Host galaxy extinction: There has been growing evidence that dark bursts are obscured in 
their host galaxies e.g. [Q— Rj. The detection of the near-infrared (IR) OAs of some GRBs 
(which would have been considered as dark bursts because their OAs were not detected in any 
bluer bands) provides evidence for dust obscuration [||]. These studies generally show that 
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GRBs originating in very red host galaxies always show some evidence of dust extinction in 
their afterglows. Also a significant fraction of dark burst hosts have extinction columns with 
Ay ~ 1 mag, and some as high as Ay = 2 — 6 mag |^]. Rossi et al. [^] performed a search for 
the host galaxies of 17 bursts with no optical afterglow. They find in seven cases extremely 
red objects in the error circles, at least four of them might be dust-enshrouded galaxies. 



2. Malmquist bias for GRB redshift measurement 

A telescope limiting magnitude (mjj defines a threshold for obtaining a redshift. Because of 
the different instruments engaged in spectroscopic redshift measurement it is difficult to constrain 
this important parameter. For defmiteness, we take = 24, approximating an average limit for 
the telescopes used to acquire the redshifts for the faintest bursts. 

With a reference time t c of 1 day, the OA limiting luminosity can be calculated at time T z [|l0|], 
assuming all OAs fade in luminosity by a °< ?~ ! : 

M L (z) =m L - 51og 10 (^) - ^log 10 (|) + 8m(z). (2.1) 

We use T z pa 400 min, which is the mean time taken to acquire a redshift [[|], and 8m(z) accounts 
for a possible increase in luminosity with redshift from either intrinsically brighter sources or less 
dust obscuration in the hosts. Using the above definition, and OA luminosity function (p(M), the 
OA selection function defines the fraction of observable OAs with redshift: 

¥oa{z) = / (p(M)dz . (2.2) 

Jm l (z) 

3. Results 



Equation ^£ can be linearly combined with other completeness functions for dust extinction 
and the redshift desert. For the latter, we employ simple scaling functions that reduce the complete- 
ness over certain redshift ranges. The GRB redshift probability distribution function, that includes 
the above selection effects, can be expressed as: 

dViz) e(z) 

P(Z)=C— — — ^V^Swift (Z)VOA(Z) ^Desert (z)Vbust(z) (3.1) 

dz (1+z) 

where C is a normalization constant, dV(z)/dz is the volume element for a flat-Lambda cosmology 
and e(z) is the dimensionless source rate density evolution function (scaled so that e(0) = 1). We 
assume that e(z) tracks the star formation rate history. 



4. Summary 

Figure [j] plots the expected redshift distribution including these completeness functions. The 
optimal model based on a Kolmogorov Smirnov (KS) probability > 95% is one that includes se- 
lection effects. Dust extinction produces a (30 — 40)% reduction in redshifts in z = — 3. The 
fraction of missed redshifts from the reshift desert is < 25% in z = 1.5 — 2.5. At high-z the OA 
Malmquist bias results in 20% of redshifts missed in z = — 5 and increases to 50% missed out 
to z = 10. 
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Figure 1: The curves plot the expected redshift distribution for models with and without optical selection 
effects, and the normalized distribution of 136 Swift triggered absorption spectroscopic redshifts from the 
OA. The model curve (with selection effects) is constrained by a KS probability > 95%. 
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